Excitons in QED 3 and spin response in a phase-fluctuating d-wave superconductor 
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We study the particle-hole exciton mode in the QED3 theory of a phase-fluctuating d-wave su- 
perconductor in ladder approximation. We derive a Schrodinger-like equation for the exciton bound 
state and determine the conditions for its existence. We find the dispersion of this mode below the 
particle-hole continuum and compare our results with the resonance observed in neutron scattering 
measurements in cuprates. 
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I. INTRODUCTION 

Strong electronic correlations in high-temperature su- 
perconductors are believed to underlie their rich phe- 
nomenology. The superconducting gap in these materials 
exhibits ci-wave symmetry with four nodes. Since the gap 
vanishes and changes sign at these nodes, the low-energy 
fermionic excitations have linear, Dirac-like, dispersion. 
An important probe of the correlations of these nodal 
quasiparticles is provided by inelastic neutron scatter- 
ing measurements, which give us information about the 
momentum and energy dependence of the correlations 
of electronic spin degrees of freedom through their cou- 
pling to scattering neutrons. Over the past decade a 
consistent experimental picture of these correlations has 
been emerging in different families of cuprates, so far 
the exclusive hosts of high-temperature superconductiv- 
ity. These include the double-layer [u |2|, p i [1 , |a, a LZ I 
Bi- and Y- and the single-layer [1, E IS till ElEl I 
La- and Tl-based families. The salient features of this 
picture are the following: (1) A resonance peak at the 
antiferromagnetic ordering wave vector, (7r,7r), at a res- 
onance energy uj res (« 41 meV at optimal doping); (2) 
A two-dimensional [14| incommensurate structure below 
the resonance energy with maxima located at (tt ± 5, tt) 
and with a possible weak inward dispersion toward (n, tt) 
as lo — > aj res ; (3) An incommensurate structure above the 
resonance with outward dispersion away from (tt, tt) for 
higher energies. The above dispersions have been taken 
to suggest that the commensurate and incommensurate 
peaks may have a common origin. 

In addition, the response is usually not discernible from 
the background below a certain energy, referred to as the 
"spin gap." However, for the sake of clarity, we reserve 
the label "spin gap" strictly for the minimum value of 
particle- hole continuum, oj sg , and call the former "spin 
response threshold." Thus, whether the observed sig- 
nals at (tt, tt) and w rcs , and/or incommensurate peaks 
are true resonances is closely linked to the question of 
whether they occur below the corresponding spin gap. 
Accordingly, one could think of two possible theoretical 
scenarios. 
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In the first scenario, there is no true resonance in the 
system, that is, no particle-hole bound state, and the 
peaks are merely the maxima of the spin response. The 
nodal quasiparticles will have a finite spin gap every- 
where except at the nodes. At (tt, tt), the peak would 
be interpreted as the overlap of the four incommensurate 
responses, one for each node, at the center. The spin 
gap and the spin response threshold are the same in this 
scenario. Usually this should also mean that the cen- 
tral peak falls off slowly as l/u>, rather than being sharp. 
However, as noted above, the central peak observed in 
experiment seems to be very sharp, in some cases limited 
only by the energy resolution of the measurement. 

In the second scenario, there would be a true resonance 
below the spin gap given by a (5-function in the spin 
response, due to the formation of particle-hole bound 
states — the so-called spin excitons. In this scenario, the 
spin gap is distinct from the spin response threshold. It 
would still be important to determine the continuum re- 
sponse so we could decide what effects derive from which 
source. However, it would also be very important to de- 
termine the existence (or lack of) and the properties of 
such excitons in a candidate theory of cuprates. 

In this paper, we study the existence of spin exci- 
tons within the QED3 effective theory of underdoped 
cuprates [HI, [H, EE EH ■ The QED 3 theory describes a 
(i-wave superconductor in which the phase of the super- 
conducting order parameter is fluctuating. These phase 
fluctuations are encoded in a pair of U (1) gauge fields, the 
Berry gauge field and the Doppler gauge field, that cou- 
ple to quasiparticle's spin and charge degrees of freedom, 
respectively. Both gauge fields are massive in the su- 
perconducting phase, which, accordingly, exhibits sharp 
nodal excitations. As the system is underdoped, however, 
the phase-fluctuations grow stronger and eventually de- 
stroy the phase coherence of the superconducting order 
parameter. The resulting state may be shown to be in- 
sulating [l9|. Whereas outside the superconductor the 
Doppler gauge field remains massive, the Berry gauge 
field becomes massless. This allows the general chiral 
symmetry breaking instability in the QED3 [20L |2"H to 
become operative, which in the present context implies 
the antiferromagnetic ordering in the system [HI [2^ . 

Herbut and Lee [H[ have previously calculated the 
QED3 spin response in the superconducting phase in 
a low-energy approximation and found no exciton res- 



onance. In their approach the Berry gauge-field propa- 
gator in the superconducting phase was approximated by 
a constant, equal to its infrared mass. This result falls 
therefore into the first scenario discussed above. This 
"constant-mass" approximation is valid below a momen- 
tum cut-off that is proportional to the gauge-field mass. 
So, it becomes less reliable in the underdoped region 
where the mass decreases and eventually vanishes at the 
underdoped transition. In the present paper we go be- 
yond the constant-mass approximation by including the 
momentum-dependence of the inverse propagator of the 
Berry gauge field, which becomes linear at high momenta. 
We do so by considering the scalar vertex for the spin re- 
sponse in a ladder approximation and by reducing the 
problem to an approximate Schrodinger equation that 
describes the formation of particle-hole bound states in 
the cJ-wave superconductor. We find that, for a strong 
enough coupling between the nodal quasiparticles medi- 
ated by the Berry gauge field such bound states exist. 
We derive their dispersion and compare our results with 
the experimental picture outlined above. However, due 
to the inherent gauge dependence of our ladder approx- 
imation we cannot definitely determine whether, for the 
physical values of parameters, QED3 theory of cuprates 
is in this strong-coupling regime. 

We note that the resonance is absent in La-cuprates. 
This could be due to different physics from competing 
stripe order that is known to exist in this family and re- 
sults in a smaller critical temperature, T c . Alternatively, 
within our theory, one might expect the resonance to be 
harder to discern from the incommensurate structure for 
a smaller T c . See Sec. |V] for a discussion. 

The paper is organized as follows: In Sec. Hi] we formu- 
late the ladder approximation for the scalar vertex and 
the spin response. In Sec. IIIII we derive the Schrodinger 
equation for excitons, find the resulting expression for 
resonant spin response, and determine the conditions for 
their existence. In Sec. IIVI we discuss the dispersion of 
the excitons and compare with experiments and existing 
literature. We summarize our results in Sec. [V] where 
we also comment on the issue of gauge dependence men- 
tioned above. The details of some of our calculations are 
given in two appendices. 



II. LADDER APPROXIMATION 

The spin part of the effective action for the nodal quasi- 
particles in the fluctuating d-wave superconductor [191 ] is 
given by 
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FIG. 1: Ladder approximation for the spin response 
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FIG. 2: Ladder approximation for the scalar vertex, T(k,p). 



The Dirac spectrum is assumed for quasiparticles near 
the nodes with anisotropic velocities t>i l(tl = (1,vf,va) 
and — (1,va,vf)- Here vf/va ~ 10 is the ratio of 
the Fermi velocity to the gap gradient at the node. We 
will set vf — va throughout our calculations and restore 
their values by rescaling the corresponding momenta at 
the end. Only the transverse components of the gauge 
field enter the action in Eq. |T]): The gauge-field propa- 
gator has the form D^{p) = (8^ - (1 - Qp»p v )D(j?), 
where El 



7T , 771 . 

F(z) = (4z 2 + l)tan" 1 — - 2z 



(2) 
(3) 



where the number of flavours, N — 2, is the number of 
pairs of nodes. The bare value of the coupling g is unity. 



and £ signifies a (nonlocal) gauge-fixing term [24]. The 
gauge field has a mass Z? -1 (0) = I2m/ir. At high ener- 
gies, D~ 1 (p) ~ \p\, which is an exact result [2^|. Most of 
what we will say about excitons in QED3 does not de- 
pend on the exact form of the interpolation between the 
infrared mass and the ultraviolet linear dependence of 
D~ 1 (p), but for the sake of definiteness we will present 
our results for this particular form of the propagator, 
which is found from the dual theory of Refs.[26] andll9l. 

The spin operator, S z , is related to the Dirac fields 
through [lfl, 

^•(r)V i (r) = 4cos(2K i T)5 2 (r) I 

where K^a = (±kF,kp) denote the positions of two 
of the the nodes in the Brillouin zone. We as- 
sume spin-rotational symmetry. Thus, x /il/ (k, u>) = 
(^(-k, -w)S v (k,w)) = Spvxfcu), where x( k , w ) = 
{S z (-k, -uj)S z (k,oj)) . Writing k = 2K, + p, we find 

x{p) = T6 / d3re ~*" r E<^w^( r )^( )^(°)>- ( 4 ) 

ij 

In order to calculate the spin response in the QED3 
action we will adopt a ladder approximation for the four- 
point correlator shown diagrammatically in FIG. [1] 
The first diagram represents the "bare" spin response in 
the absence of gauge-field interactions. We denote it by 
XoO)- 

It is standard to reformulate the ladder approximation 
in terms of gauge-field-induced corrections to the scalar 
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vertex for the interaction between the external source 
(neutrons) and spinons. These vertex corrections are 
shown in FIG. [2] Defining the scalar vertex, T(k,p), as 
the amputated diagram (without the external legs) on 
the left-hand side, we have 



N f 

x(p) = - TE J 



d 3 k 
(2^)3 



ti[Go(k)T(k,p)Go(k+p)], (5) 



where Go(k) = —ij^k^/k 2 is the bare four-component 
spinon propagator in the superconducting state. From 
FIG. [2|we see that the scalar vertex satisfies the following 
Bethe-Salpeter equation, 



T(k,p) = t-g 2 



d 3 g 
(2tt) 3 



j fl G {q)r(q,p)G (q + p) 

x^D^ik ~ q). (6) 



We may use the symmetries of the QED3 action to 
expand the scalar vertex in terms of form factors [27| , 



T(k,p) = tF 1 (k,p) + (r flv k fl p v F 2 (k,p) 

^l^k^F^k^) + l^P^F^k^p). 

Here, a^y = 5[7 M ,7i/]. Then, Eq. ^ is reduced to 
d 3 q 



(7) 



Fi(k,p) = 5a + 



A 



(2tt) 3 ^ 



K iJ (k,q 1 p)F J (q,p), (8) 



with the matrix ify having a block-diagonal form in 
terms of two two-by-two matrices. For the reasons ex- 
plained below the most important element is 



Kn(k,q,p) 



(2 + 09 2 q 2( iq ^i D(k- q ). (9) 
q {q +p) 



The block-diagonal form of implies that the form 
factors F 3 and F4 completely decouple from Fx and F 2 . 
Together with Eq. (6), which implies that F3 and F4 do 
not contribute to the function this allows us then to 
neglect them altogether for the purpose of calculating the 
spin response. 

Within the ladder approximation for the scalar vertex, 
our manipulations have so far been exact. At this point, 
we notice that the dominant contribution to x comes 
from F%, which is of order g . Since F 2 ~ g 2 , we will 
then neglect F 2 to obtain the set of equations [28l l29j , 



f ^ N 
x(p) = -j 



d 3 k k-(k+p) 



(2tt) 3 k 2 (k+p) 
f d 3 , 

Fl(fc ' p) ^ 1 + A i(^)3 n 9+P y 



2 Fi(k,p), 
d 3 q q-(q+p) 



(10) 

D(k-q)F 1 (q,p), 
(11) 



for the spin response, with A = (2 + £)g 2 as the (gauge- 
dependent) coupling strength. These equations provide 
the basis for our further study of the exciton modes in 
QED 3 . For a constant D(k—q) they reproduce the results 
of Ref . HI 



Note that Eqs. (fTTj) and (fTO)) explicitly depend on the 
choice of gauge. This arises because the ladder diagrams 
in FIG. O are not gauge-invariant if the bare fermion 
propagator is used. The gauge invariance can be re- 
stored by choosing the (non-local) gauge-fixing param- 
eter ^ so that the bare fermion propagator satisfies the 
usual Ward-Takahashi identity with the vertex. In prin- 
ciple this procedure will yield a momentum-dependent 
function, £ (p) (29j . With m — 0, we find Nash's gauge 
£ = 2/3 [3Cj] whereas in the opposite limit of constant- 
mass approximation one finds Feynman's gauge, £ = 
1. In the rest of the paper, we choose to work with 
a momentum-independent gauge-fixing parameter £ for 
simplicity. 



III. EXCITONS IN QED a 

Our strategy to solve the vertex function is to obtain 
a (set of) differential equation(s) from the integral equa- 
tion (fTTj) . This is similar to the classic derivation of exci- 
ton bound states in metals and semiconductors [3lL |32| . 
We will show that there is an approximate Schrodinger- 
like equation that captures the bound-state content of 
the vertex function. Similar Schrodinger equations have 
been derived before for the non-relativstic limit of vertex 
functions in QED [H, Based on this Schrodinger 

equation, we will determine the conditions for the exis- 
tence of the exciton bound states. 



A. Schrodinger equation 



Let us define a generalized response 

d 3 fc k^{k+p} v ik „ 



^ f d ' 



(2tt) 3 k 2 (k+p) 2 



Fi(k,p), (12) 



where ay, = (o\, a 2 , 03) are Pauli matrices (up to cyclic 
permutations). Then Eqs. (fT0j) and (fTTj) may be written 
as 

N N 
X (p) = T tr[$(O,p)] = T 0(O,p), (13) 

Fi(q,p) = 1 + \J d 3 re- iq - r D{r)(j){r 1 p). (14) 

Note that r is a real-space variable and D(r) indicates 
the real-space gauge-field propagator. 
We see that 

-[(ip+pxcr)-d+d-d)]<f>(r,p) = 5 {3) (r)+\D(r)tr [$(r,p)] . 

(15) 

We may shift away the momentum p (which enters the 
equation only as a parameter) in the derivative by com- 
pleting the square. This is achieved by a phase transfor- 
mation of the form, 



$ — > $' = exp 



^(ip+px a) -r 



(16) 
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We note that <f>'{0,p) 
tion is, then, 



<f)(0,p). The transformed equa- 



[~d-d+(p/2) 2 ]^' -\D{r)e^ L ti[ 



(17) 

Where we have denoted the orbital angular momentum 
of the particle-hole system by L — r x p. This is an exact 
mapping of the integral equation (TTTJ to a set of coupled 
Schrodinger-like differential equations. 

Even though we have succeeded in deriving a set of dif- 
ferential equations for our vertex function, they are still 
very complicated. However, for the purpose of study- 
ing the bound-state content of the vertex function, we 
may simplify these equations further down to a single 
Schrodinger-like equation. A numerical solution seems to 
be the only useful alternative. As shown in Appendix YK[ 
this approximate, decoupled equation is given by 

[(p/2f - d ■d-V eB (r,p,m)]<f>'(r,p) = *< 3 >(r), (18) 

where, for < \p\ < 2m, 

Asin 2 6» s 



V bB = A 1 



16 



cosh^ \L/2\ D(r). 



(19) 



Here, 8 = cos _1 (p ■ f). We expect Eq. (fTS)) to provide a 
valid description of the bound states for an adequately 
short-ranged interaction D(r). This is the case when 
< \p\ < 2m. For \p\ > 2m, Eq. §T§§ is valid only at 
short distances. At large distances r ^> (\p\ — 2m)" 1 , the 
potential becomes infinite for 6 > sin _1 (2m/|p|): 



A 2 

o 

p 



sinh 2 \L\ D 2 (r) -> -co. 



(20) 



We note that in the case of constant-mass approxima- 
tion [23[ for which D{r) ~ S^(r), and also for p = in 
the general case, Eq. (fl~8|) follows exactly from Eq. (fT7|) 
with Kff(r) = XD(r). 



B. Resonant spin response 

We may now find x{p) by solving the following three- 
dimensional Schrodinger equation for the eigenvalue 
e„(|p|, m) and the normalized cigenfunction ip n {r, to), 



[-8 • d - V e s] Ipn = e n ^n- 



Then, 



X(p) = x^(°'f) =iV E 



\M0,\p\,m)\ 2 
4e n (\p\,m) +p 2 



(21) 



(22) 



By analytically continuing to real frequencies, pa — > 
—iu> + + , and denoting pu = \/p 2 — , we find that 
since Se„(pM,"i) = the imaginary part of the spin 
response observed in experiment is given by 



where the sign ± = sgn(o;)sgn (l + de n (pM, m)/dp 2 yI ). 
If e n (pM, Tn) were not real there would be no (5-function, 
hence no resonance. 

From Eq. (|23[) . we see that the necessary and sufficient 
condition for the existence of excitonic resonances at p 
and uj < |p| is that our Schrodinger equation (f2T|) admits 
bound state solutions with real and negative eigenvalues, 
that solve the equation 



( \ Pc 



(24) 



We will now study the existence of such bound states. 



C. Existence of excitons 

As shown in Appendix [B] the gauge-field propagator 
D(r) scales as 1/r 2 . For \p\ > 0, we may rescale the space 
as \p\r = z to see that the energy spectrum satisfies the 
scaling relation 



e&(|p|,«i) =p 2 e b (2m/\p\). 



(25) 



The scaling function e b is the bound-state energy eigen- 
value for the rescaled Schrodinger equation, 



V cS (z,p) ipb{z,[J.) = e b {n)tpb{z, p), (26) 



where \i = 2m/\p\ and the rescaled potential is 



(27) 



The potential V e g has an inverse-square form for small z. 
The condition (|2~4"|) for the resonant response is satisfied 
for e b (p c ) = -1/4; then, p c = 2m/fi c . 

The inverse-square potential is an instance of confor- 
mal anomaly, i.e., the breakdown of scale symmetry in 
quantum mechanics. The potential — A/16z 2 has an in- 
finite number of negative energies for A > Aq = 4 and 



is unbounded from below [35|. Thus, in this "strong- 
coupling" regime, the problem needs to be renormalized 
for it to be physically meaningful. Although there are 
different ways of doing so, the result is unique [36]. The 
inverse-square potential has a single renormalized bound- 
state. The energy of the renormalized bound-state can- 
not be determined within the theory; it is an input of 
the theory either from experiment, or from the physics 
at higher energies, beyond the domain of physical validity 
of the potential. The continuum spectrum that is renor- 
malized into a single bound-state is called the "conformal 
tower." The bound-state wave function is given by 



il>o{z) 



k 3 K {k\z\) 



M ) 



Nj2±\M0,PM,m)f 



(pIi 



^n(j>M TO ))(23) wn ere — k is the renormalized bound-state energy, and 
Kq is the modified Bessel function of the second kind. 



For \x > 1 the potential V e s decays exponentially at 
large distances. However, the inverse-square form at 
small z is expected to be sufficient for the existence 
of a conformal tower in the strong coupling A > A*. 
By using ?/>o as a trial wave function to calculate the 
energy one can check that this expectation is met for 
A*/Aq 1.268. The corresponding critical charge is given 
by g c w 2.252/ v / 2+~£. In Nash's and Feynman's gauge, 
for example, one finds g Cy N = 1-379 and g c .p — 1.300, 
respectively. The strong-coupling regime is found for 
9 > 9c- 



For < \i < 1, the potential V c g 



^2(sin9-iJ.)\z\ 



—oo for 9 > sin /i. So it is less clear whether the con- 
formal tower would exist in this limit. However, the same 
way as above and by choosing n > (1 — /z)/2, we have 
checked that the conformal tower still persists for the 
same strong-coupling regime. This result reflects the fact 
that the conformal tower for bound-states is essentially 
a short-distance phenomenon produced by the singular 
behaviour of the potential at the centre. 



IV. DISCUSSION AND COMPARISON WITH 
EXPERIMENT 

In this section we will assume that QED 3 theory is 
in the strong-coupling regime characterized in Sec. IIII CI 
and investigate the consequences of this assumption for 
the spin response observed in neutron scattering mea- 
surements. 

The dispersion of the exciton mode can be read off 
from Eq. {23]) and Eq. {24]) to be 



*(p) 



; (P) -fc 



(28) 



where w sg (p) is the "spin gap," given by the minimum 
of the particle-hole continuum. With the linear Dirac 
spectrum for spinons, it is 



^sg(p) 



2 2 
VpPx 



2 2 



where we have now restored the Fermi and gap veloci- 
ties. In this formula, p x \ \vf is the component of the wave 
vector measured from the given node in the nodal direc- 
tion and p a ||wA is the one in the perpendicular direction. 
There is one dispersion branch for each node. Due to 
the strong anisotropy va/vf "C 1, the spin gap is much 
more sensitive to the changes of p in the nodal direction. 
This causes the dispersion (|2"5]) to assume a characteristic 
shape shown in FIG. [3] This is in qualitative agreement 
with experiment if we identify the exciton mode with the 
observed resonance peak [l(J. We will now discuss some 
features of the dispersion. 

In the diagonal direction, FIG.[3Ja), there are three ex- 
citon branches. Each of the nodes on the same diagonal 
line contribute one branch, and the other two nodes pro- 
duce overlapping branches. They all cross at (jr, it) and 
hence a strong resonance is expected here. The overlap- 
ping branch is rather flat. Its curvature is determined by 










ik 


(b) v 









-0.20 



-0.11 



0.11 



0.20 



FIG. 3: (Color online) The dispersion branches of the exciton 
in (a) the diagonal direction; and (b) the parallel direction, 
as shown in the insets. The wave vectors are in reciprocal 
lattice units, and the origin is at (n,n). The value of param- 
eters used are va/vf = 0.15, fci? = OAn and p c /vF = tt/16. 
The continuum spin gap is shown by the thick blue line and 
vanishes at the nodes. 



the ratio va/vf- It quickly enters the continuum away 
from (7r,7r) and is, then, presumably damped. The other 
two branches fall off rather steeply (due to the square 
root in w res ) and terminate a distance p c /vf away from 
the corresponding nodes. 

In the parallel direction, FIG.[3](b), there are only two 
sets of overlapping branches, each from a pair of nodes 
that map to each other upon reflection about the paral- 
lel line. Here, too, we observe a crossing at (tt, it) and, 
therefore, an enhanced response. The absence of the flat 
branch means that the momentum width of the resonance 
at (-7T, 7r) is only bound by the momentum resolution of 
the experiment in this direction. Again, the modes dis- 
perse down to zero energy at an incommensurate position 
(ir ± 5 lnc , it) where 



2(tt- 2k F ) - V2p c /v F 
l + v 2 Jv 2 F 



< 2(tt - 2k F ) - V2p c /v F , 
(29) 

and the approximation is made for va/vf <C 1. 

In our Schrodinger equation approximation, the 
strength of the peak is given by |^ n (0,p c )| 2 at the reso- 
nance (pm = Pc) and is the same for different momenta. 
That is, the whole branch below the continuum spin gap 
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has the same intensity. This does not agree with exper- 
iments in which the resonance peak appears to go away 
for low energies. 

We also briefly note the effects of doping, x, on the 
dispersion of the exciton. In our formulation, doping 
enters through the dependence of the gauge-field mass, 
m(x), the position of the nodes, kp(x), and the spin gap , 
Lu sg (p,x). The gauge-field mass scales with T c [IjJ, 1 23l | . 
Thus, it decreases with underdoping and vanishes at the 
underdoped superconducting transition. As we approach 
half-filling, the nodal points move towards (tt/2, 7r/2) and 
Uf increases upon underdoping. This causes the spin gap 
to decrease [371 ] . From our scaling relations in Sec. IIII Cl 
we see that the value of p c (x) ~ m(x) also decreases with 
underdoping. Based on these trends, we can expect that 
with underdoping, and hence with decreasing T c : (i) the 
resonance energy w ros decreases and merges with the de- 
creasing spin gap; and (ii) the incommensurability 8- mc 
decreases and merges with the nodal points. These are 
in qualitative agreement with experiment [ToL l38l . l39l . [40| . 
We note, in contrast, that if kp remains fixed while 
T c changes in a set of experiments, Eq. (f25)) predicts 
that the incommensurability S inc will increase linearly 
with decreasing T c . Underdoping, however, increases kp, 
hence the available area for the incommensurate struc- 
ture shrinks. 

The resonance and the incommensurate features go 
away at the transition, where m = T c = 0. As T c — > 
the resonance merges in with the continuum response. In 
our theory, we interpret the absence of a resonance in La- 
based cuprates with an anomalously small T c compared 
to other cuprate families, as the difficulty to discern such 
a merging feature. If this is indeed what is happening, 
the same behavior should be seen in other families as 
T c — > 0. Of course, it is possible that the absence of the 
resonance in La-cuprates is due to the different physics, 
that would give rise to stripes, for example. 

To conclude our discussion, we briefly compare our re- 
sults to existing theoretical literature. The spin response 
in cuprates has been studied extensivel y El |42[ |43| . 0, 

EE Efl EE EE H M M HI HI HI Wm- Th e most 

common approach to the problem is the random-phase 
approximation (RPA) derived or postulated in the con- 
text of a spin-fermion model such as the (extended) t-J 
model. With a d-wave gap and by appropriately choos- 
ing the hopping amplitudes (fermiology) the resonance 
arises as the pole of the RPA response at a step (Van 
Hove) singularity of the bare spin response. The disper- 
sion of the resonance is the result of the dynamical nest- 
ing of the energy spectrum [50(. It has also been argued 
that the d-wave gap as opposed to the fermiology is the 
root of the resonance structure [4(| |49[ . This resonance 
has been interpreted as a spin collective mode or exciton. 
Various phenomenological parameters such as higher har- 
monics of the sup erconducting gap [5l|, interlayer cou- 
plings [HI, [H, HJ , and Fermi surface geometry [48[ have 
been modeled to fit the experimental results. 

Clearly, our work shares with the above studies the 



explanation of the resonance as a spin-collective mode. 
We calculate the spin response in the ladder approx- 
imation which is also the basis of the RPA response 
above. However, the mechanism that gives rise to the 
exciton differs in our approach. In fact there is no sim- 
ple RPA form for our response, and the bare response 
QVo(P- w ) ~ \/ <jj2 — p 2 Q(ui 2 — p 2 ) does not contain any 
step singularity [23[. The exciton we find is due to the 
strong coupling between the spinons and vortex defects of 
the superconducting order parameter through the Berry 
gauge field. Perhaps more importantly, the dispersion 
found here (FIG. [3]) persists, at least in our simplified 
analysis, to very low energies. The dispersion in the 
RPA approach has a similar shape to ours but practi- 
cally merges into the continuum at a finite frequency be- 
low the commensurate uj rcs . In most RPA studies the 
form of the spectrum away from the nodes (e.g. the cur- 
vature of the Fermi surface) is an important part of the 
explanation. We have, in contrast, only kept the lin- 
ear part of the spectrum, which is justified within our 
low-energy theory, and assumed the mass of the Berry 
gauge field is small enough for the high-energy tail of 
the gauge-field propagator to be important. Finally, the 
RPA studies are phenomenological and rely on numerical 
calculations. They reproduce rather well many aspects 
of the experiments. Our work is, on the other hand, an 
analytical study within the effective QED3 theory of un- 
derdoped cuprates. To be able to proceed analytically we 
have relied on various approximations. Although we have 
argued they are reasonable, these approximations might 
be too restrictive to contain all the physics revealed in 
neutron scattering measurements. 



V. CONCLUSION 

We have studied the problem of particle-hole bound 
states of spinons (excitons) in the QED3 effective theory 
of a phase-fluctuating d-wave superconductor. This the- 
ory contains a massive gauge field with an exact confor- 
mal propagator at high-energies. We employed a ladder 
approximation and derived an approximate Schrodinger- 
like equation for the bound states. We discussed the con- 
ditions for the existence of excitons and concluded that 
they would exist in the strong-coupling regime. We de- 
duced the dispersion of the excitons and compared our re- 
sults with neutron scattering measurements in cuprates. 

This work complements the earlier work of Herbut 
and Lee who discuss the continuum spin response 
in QED3 in a low-energy approximation. It was found 
there that dispersing incommensurate and commensu- 
rate peaks similar to the resonance structure exist in 
the continuum. The qualitative behaviour of this con- 
tinuum spectrum remains the same when higher-energy 
effects are included, but the numerical values of spin gap 
are reduced and are in fact close to the experimentally- 
observed range of resonance energies [13] ■ Thus, the con- 
tinuum spectrum must be accounted for in extracting 
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information on the exciton spectrum from experiments 
with a finite energy and momentum resolution. 

Our vertex function is by definition gauge-dependent. 
In this work we used a non-local gauge-fixing term, 
^D(p)pnp v , in the gauge-field propagator with a 
momentum-independent parameter £. However, this is 
not enough to ensure the gauge invariance of the response 
for D{p) given by Eq. @: enforcing the Ward-Takahashi 
identity for the bare spinon propagator yields [2!| a 
momentum-dependent parameter 2/3 < £(p) < 1, with 
the lower and higher bounds corresponding to the ultravi- 
olet and infrared limits, respectively. A renormalization- 
group analysis on the original dual theory of Ref.[l7|, from 
which one can derive the action in Eq. ([1]), has been per- 
formed in Ref. HE It was found that the charge g cannot 
be renormalized due to the non-analyticity of the spinon 
polarization ~ \p\. So, within our ladder approximation 
and for the bare value g = 1, it is likely that the gauge- 
invariant response is not in the strong-coupling regime. 
Further study of the gauge- invariant response seems nec- 
essary. 

More comprehensive numerical studies of this prob- 
lem will be very useful. They can be applied, at var- 



ious stages, to the original vertex equations, or to the 
Schrodinger equation for the bound states. Finally, it 
is interesting to see whether other spin-collective modes, 
say triplets, can be formed in QED3. 
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APPENDIX A: DERIVATION OF EQ. (HHJ) 

We decompose the generalized response into symmetric 
and anti-symmetric components as $' = 10' + icr ■ <p' A . 
Then, keeping in mind that L — r x p, from Eq. (|17[) we 
find 



d-d- Acosh 2 |L/2|L>(r) 



(j)'(r,p)+i-smh\L\ D(r) L ■ cf>' A [r,p) = 5 {3) (r), 



d ■ d + A sinh 2 \L/2\D{r)LL 



(f>' A (r,p) +i-smh\L\ D(r) L <p'(r,p) = 0. 



(Al) 
(A2) 



r 



This is a set of coupled Schrodinger equations. However, 
we notice that the second equation contains a repulsive 
potential and consequently has presumably no negative- 
energy bound-states. Thus, for |p| > 0, its Green's func- 
tion drops exponentially with increasing distance. In or- 
der to account for its effects on the bound states of the 
first equation, we may then solve the second equation for 
finite p simply as 

2A 

cf> A (r,p)tt -i^rsinh.\L\D(r)L<j)'(r,p). 

By plugging this expression into the first equation we 
have [(p/2) 2 - d ■ d - V cS {r)]<j)' (r,p) = 5&(r), where 



Ves(r) = A 



4A 9 
1 - -^sinh 2 \L/2\ D(r) 



cosh 2 \L/2\D{r) 



For < \p\ < 2m, due to the exponential fall-off of 
D(r) (see Appendix IB"]), we have 



P 



■sinh 2 \L/2\ D 2 {r) 



' r X f D(r). 

lQ r 2p2 



Thus, we obtain Eq. (fig)) . For |p| > 2m we still find 
the same behaviour as above at short distances. But. for 



large distance, the ~ A 2 term overwhelms the other term 
and we have 



V eff (r > 



1 



A 2 



|p| - 2m' 
as claimed in Eq. (|20| 



sinh 2 \L\ D 2 {r), 



APPENDIX B: REAL-SPACE GAUGE-FIELD 
PROPAGATOR 

Here we derive the real-space interaction D(r). Since 
there is a gap in the low energy limit, we expect an ex- 
ponential decay in the long-distance behaviour of D{r). 
Also, the tail in the high-energy limit should trans- 
late to a 1/r 2 singularity at short distances. We start by 
noting that 



. . m 

D(r) = 

4-7rr 



A/ '2m 



~ , „ m 



F(l/2z)sin(2m|r|z)dz (Bl) 
(B2) 
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where A is an ultraviolet cutoff and 



D(r) 



8tt r 



-A/2m 



1 + z 2 



tan 1 z sin(2m|r|z)d2:. 



-A/2r. 



(B3) 

We compute D(r) using contour integration. To this 
end, we need to choose two branch-cuts to define tan -1 z. 
We take them to be [+i, -Moo) and [— i,— ioo) such that 
z-i = ne iei with < Q x < | an d z + i = r2e i9 2 w j t h 

-f < 2 < Then, 

tan -1 z = -— + log — + i(9 2 - 9i). 
2 n 

We take a contour C that includes the real axis and closes 
on itself in the upper-half plane, except that it avoids the 
upper branch-cut by tracing the path C\ = — F(0 + ) U 
F(0~) where F(s) — [e + i,e + ioo). So, we can now 
write 



DM = ZL 3 
v ' 8vrr 



m 
8r 



c Jc 



l + z A 



■ tan 



ze 2lm ^ z dz 



(B4) 



So, we find 
D(r) 



m [°° u 2 - 1 



8r u 2 



-2m|r | 



16r 2 8|r| 



—E 2 (2m\r\). (B6) 



Note that D(r) > 0. Here, £ n (a;) = e- ux du/u n is 
the exponential integral. For n — 2 it has the following 
asymptotic behaviour 

Eo(x) - J e_1 (i-? + •••). *-»°°> 

21 J \ l + ( 7B -l + loga;)a; + --- , a; 0, 

where 7b is the Euler-Mascheroni constant. It leads in 
turn to the following asymptotic behaviour for D(r): 



D(r — » oo) 
D(r -> 0) 



g — 2m | r | 

64m|r| 3 



+ 0(|r|- 4 e^ m l r l), (B7) 



1 m 



16r 2 41 



+ (9(m 2 log(m|r|)).(B8) 
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